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Abstract 


In this paper we investigate the limit behavior of the solution to quasi-static Biot’s 
equations in thin poroelastic flexural shells as the thickness of the shell tends to zero and 
extend the results obtained for the poroelastic plate by Marciniak-Czochra and Mikelic 
in [16]. We choose Terzaghi’s time corresponding to the shell thickness and obtain the 
strong convergence of the three-dimensional solid displacement, fluid pressure and total 
poroelastic stress to the solution of the new class of shell equations. 

The derived bending equation is coupled with the pressure equation and it contains 
the bending moment due to the variation in pore pressure across the shell thickness. The 
effective pressure equation is parabolic only in the normal direction. As additional terms 
it contains the time derivative of the middle-surface flexural strain. 

Derivation of the model presents an extension of the results on the derivation of 
classical linear elastic shells by Ciarlet and collaborators to the poroelastic shells case. 

The new technical points include determination of the 2x2 strain matrix, independent of 
the vertical direction, in the limit of the rescaled strains and identification of the pressure 
equation. This term is not necessary to be determined in order to derive the classical 
flexural shell model. 

Keywords: Thin poroelastic shell, Biot’s quasi-static equations, bending-flow 
coupling, higher order degenerate elliptic-parabolic systems, asymptotic methods. 

AMS classcode MSC 35B25, MSC 74F10, MSC 74K25, MSC 74Q15, 

MSC 76S. 

1 Introduction 

A shell is a three dimensional body, defined by its middle surface S and a neighborhood of 
a small dimension (the thickness) along the normals to it. The shell is said to be thin when 
the thickness is much smaller then the minimum of its two radii of the curvature and of the 
characteristic length of the middle surface S. 

The basic engineering theory for the bending of thin shells is known as Kirchoff-Love 
theory or Love’s first approximation. The equations were derived by the so called ’’direct 
method” (see [T8| and references therein) and not from the three dimensional equations. A 
derivation from the three dimensional, at the rigor of the continuum mechanics, is due to 
Novozhilov and we refer again to [l 8 ] for both linear and nonlinear models. 

A different approach to deriving the shell equations is to suppose that the middle surface 
S is given as 5 = X(Zi;), where w C be an open bounded and simply connected set with 
Lipschitz-continuous boundary duj and X : tJ —)■ is a smooth injective immersion (that is 
X G and 3x2 matrix VX is of rank two). The vectors Sia{y) = '9oX(?/), a = 1,2, are 
linearly independent for all y G ZJ and form the covariant basis of the tangent plane to the 
2-surface S. 

Then the reference configuration of the shell is of the form r(n ), e > 0, where = 
uj X (—e/2, e/2) and 

r = r(y,X 3 ) = X(y) -bX 3 a 3 (y), a 3 (y) = ^ 

|ai(y) X a 2 (y)| 

The associated equations of linearized three dimensional elasticity are then written in curvi¬ 
linear coordinates with respect to {yi,y 2 ,X 3 ) G 12 . Their solutions represent the covariant 
components of the displacement held in the reference conhguration r (12 ). 

Then Ciarlet and collaborators developed the asymptotic analysis approach where the 
normal direction variable X 3 G (—e/2, e/2) was scaled by setting 2/3 = x^/e. This change of 
variables transforms the PDE to a singular perturbation problem in curvilinear coordinates 
on a hxed cylindrical domain. With such approach Ciarlet and collaborators have established 
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the norm closeness between the solution of the original three dimensional elasticity equations 
and the Kirchhoff-Love two dimensional flexural and membrane shell equations, in the limit 
as e —>■ 0. For details, we refer to the articles [S] and m and to the books and [7]. For 
the complete asymptotic expansion we refer to the review paper m by Dauge et ah Further 
generalizations to nonlinear shells exist and were obtained using F— convergence. We limit 
our discussion to the linear shells. 

In the everyday life we frequently meet shells (and other low dimensional bodies) which are 
saturated by a fluid. Many living tissues are fluid-saturated thin bodies like bones, bladders, 
arteries and diaphragms and they are interpreted as poroelastic plates or shells. Furthermore, 
industrial filters are an example of poroelastic plates and shells. 

Our goal is to extend the above mentioned theory to the poroelastic shells. In the case 
of the poroelastic plates, derivation of the mathematical model was undertaken in |16] . As 
in the case of plates, these are the shells consisting of an elastic skeleton (the solid phase) 
and pores saturated by a viscous fluid (the fluid phase). Interaction between the two phases 
leads to an overall or effective behavior described by the poroelasticity equations instead of 
the Navier elasticity equations, coupled with the mass conservation equation for the pressure 
held. The equations form Biot’s poroelasticity PDFs and can be found in m, la and in the 
selection of Biot’s publications [23] . 

The effective linear Biot’s model corresponds to the homogenization of the complicated 
pore level huid-structure interaction problem based on the continuum mechanics hrst prin¬ 
ciples, i.e. the Navier equations for the solid structure and by the Navier-Stokes equations 
for the how. Small deformations are supposed and the interface between phases is linearized. 
The small parameter of the problem is the ratio between characteristic pore size and the 
domain size. If, in addition, we consider a periodic porous medium with connected huid and 
solid phases, then the two-scale poroelasticity equations can be obtained using formal two- 
scale expansions in the small parameter, applied to the pore level huid-structure equations. 
For details we refer to the book m, the review [T] and the references therein. 

Convergence of the homogenization process for a given frequency was obtained in [19] . 
using the two-scale convergence technique. Convergence in space and time variables was 
proved by Mikelic and collaborators in m and M- The upscaling result was presented 
in detail in [16] and we avoid to repeat it here. We point out that the upscaled model 
depends on the particular time scale known as Terzaghi’s time Tt- It is equal to the ratio 
between the viscosity times the characteristic length squared and the shear modulus of the 
solid structure multiplied by the permeability. If the characteristic time is much longer than 
Tt than flow dominates vibrations and the acceleration and memory effects can be neglected. 
The model is then the quasi-static Biot model. For its derivation from the first principles 
using homogenization techniques see m- 

For the direct continuum mechanics approach to Biot’s equations, we refer to the mono¬ 
graph of Coussy m- Using a direct approach, a model for a spherical poroelastic shell is 
proposed in [22]. 

In this paper we follow the approach of Ciarlet, Lods and Miara, as presented in the 
textbooks [6] and [7], and rigorously develop equations for a poroelastic flexural shell. 

Successful recent approaches to the derivation of linear and nonlinear shell models use 
the elastic energy functional. In our situation, presence of the flow makes the problem quasi¬ 
static, that is time-dependent, and non-symmetric. The equations for the effective solid 
skeleton displacement contain the pressure gradient and have the structure of a generalized 
Stokes system, with the velocity held replaced by displacement. Mass conservation equation 
is parabolic in the pressure and contains the time derivative of the volumetric strain. 

We recall that the quasi-static Biot system is well-posed only if there is a relationship 
between Biot parameters multiplying the pressure gradient in the displacement system and 
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the time derivative of the divergence of the displacement in the pressure equation. We were 
able to obtain in m the corresponding energy estimate. Similar estimates, for the equations 
in the curvilinear coordinates, will be obtained here. 

In addition, there exists a major difference, with respect to the limit of the normalized 
633 term, compared to a classical derivation of the Kirchhoff-Love shell. In our poroelastic 
case, the limit also contains the pressure field. Furthermore, the pressure oscillations persist 
and we prove the regularity and uniqueness for the limit problem. As in the poroelastic plate 
case, it has a richer structure than the classical bending equation. We expect more complex 
time behavior in this model. 


2 Geometry of Shells and Setting of the Problem 

We study the deformation and the flow in a poroelastic shell = r(0£), L,£ > 0, where the 
injective mapping r is given by (jl.ip . for X3 G (—.£/2,.£/2) and {yi,y2) G wl, diam (lol) = L. 
We recall the middle surface S = is the image by a smooth injective immersion X of 

an open bounded and simply connected set ul C M^, with Lipschitz-continuous boundary 
dujL- We use the linearly independent vectors aa{y) = da^{y), a = 1,2, to form a covariant 
basis of the tangent plane to the 2-surface S. The contravariant basis of the same plane is 
given by the vectors a°‘{y) defined by a"(?/) • a,g(y) = 6^. We extend these bases to the basis 
of the whole space by the vector a 3 given in dni) (a3 = a 3 ). Now we collect the local 
contravariant and covariant bases into the matrix functions 


Q = [ a^ a^ a^ ], Q ^ 



( 2 . 1 ) 


The hrst fundamental form of the surface S, or the metric tensor, in covariant Ac = (oq/j) or 
contravariant = (a"^) components are given respectively by 

Oo/? = a^ • a,g, = a“ • a^, a,f3 = l, 2 . 


Note here that because of continuity of A'^ and compactness of ujl, there are constants 
> mP > 0 such that 


rrfyi ■ x < A'^(y)x • x < M'^x • x, x G y G TDl- 


( 2 . 2 ) 


These estimates, with different constants, hold for Ac as well, as it is the inverse of A'^. The 
second fundamental form of the surface S, also known as the curvature tensor, in covariant 
Be = (^a/?) or mixed components B = (ha) are given respectively by 

2 

Kfi = a^ • djiaa = ■ a«, b^ = '^ a^'^b^a, a, ^ = 1, 2. 

K,= l 

The Christoffel symbols are defined by 


T);^ = a"'• 9^aa =- 9 / 3 a'" • a«, a,/3,K; = l,2. 

We will sometime use Tj)^^ for bay- The area element along S is '/ady, where a := det Ac. 
By (12.2p it is uniformly positive, i.e., there is rua > 0 such that 

0 < ma < a(y), y G ul - (2.3) 
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We also need the covariant derivatives which are defined by 

2 

a,P,K = 1,2. (2.4) 

r=l 


In order to describe our results we also need the following differential operators: 

I 2 

7a/3(v) = -{daVfi + dpVa) - Y o, /3 = 1, 2, (2.5) 

K,= \ 

2 2 2 2 2 

Pa/3(v) = - Y + Y - Y “ Y 

K.= l K.= l r=l K.= l T=1 

2 2 


+ a,/? = 1 , 2 , ( 2 . 6 ) 

h;= 1 ft=l 

2 2 

^a/3|/3 = dp'^aP + ^ ^ 4" ^ ^ /4 = Ij 2, 

K =1 K =1 

2 

la/3 — (^a/3 1 /?) 4“ ^ ^ rQ,^7lQ,^|^ , O, /3 — 1 ) 2 , 

K=1 


defined for smooth vector helds v and tensor fields n. 

The upper face (respectively lower face) of the shell is = r{ujL X {x 3 = il2\) = 
r(S£) (respectively = r(a;L x {xs = —tjiy) = r(S)(^). is the lateral boundary, 
r£ = r(cla;i;, x = r(r|^). We recall that the small parameter is the ratio between 

the shell thickness and the characteristic horizontal length is e = ^/L <C 1. 


Table 1: Parameter and unknowns description 


SYMBOL 

QUANTITY 


shear modulus (Lame’s second parameter) 

A 

Lame’s first parameter 

Pg 

inverse of Biot’s modulus 

a 

effective stress coefficient 

k 

permeability 

V 

viscosity 

L and i 

midsurface length and shell width, respectively 

e = t/L 

small parameter 

T = r]L'i/{k^i) 

characteristic Terzaghi’s time 

U 

characteristic displacement 

P = Ufi/L 

characteristic fluid pressure 

U = {UI,U2,U3) 

solid phase displacement 

P 

pressure 


We note that Biot’s diphasic equations describe behavior of the system at so called Terzaghi’s 
time scale T = rjL'^/{k^), where Lc is the characteristic domain size, rj is dynamic viscosity, 
k is permeability and ^ is the shear modulus. For the list of all parameters see Tabled] 
Similarly as in m, we chose as the characteristic length Lc = i, which leads to the 
Taber-Terzaghi transversal time T^ab = r]i'^/{k^). Another possibility was to choose the 
longitudinal time scaling with Ti^^g = ^T^/(/c/i). It would lead to different scaling in (|2.9p 
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and the dimensionless permeability coefficient in (13.31) would not be but 1. In the context 
of thermoelasticity, one has the same equations and Blanchard and Francfort rigorously 
derived in [S] the corresponding thermoelastic plate equations. We note that considering the 
longitudinal time scale yields the effective model where the pressure (i.e. the temperature in 
thermoelasticity) is decoupled from the flexion. 

Then the quasi-static Biot equations for the poroelastic body take the following di¬ 
mensional form: 


a = 2/ie(u) -|- (A div u — ap)l in 
— div a = — /r A u — (A -|- ^) V div u-|-aVP = 0 in 

d k ~ 

-^WgP + a div u) - Ap = 0 in 

at rj 


(2.7) 

( 2 . 8 ) 

(2.9) 


Note that e(u) = sym vu and a is the stress tensor. All other quantities are defined in Table 

^ ~ k dp 

We impose a given contact force ai/ = and a given normal flux-—— = Vl at 

rj 0 x 5 

X 3 = ±^/ 2 . At the lateral boundary T^ we impose a zero displacement and a zero normal 
flux. Here 1 / is the outer unit normal at the boundary. At initial time t = 0 we prescribe the 
initial pressure 

Our goal is to extend the Kirchhoff-Love shell justification by Ciarlet, Lods et al and by 
Dauge et al to the poroelastic case. 


We announce briefly the differential equations of the flexural poroelastic shell in dimen¬ 
sional form. Note that our mathematical result will be in the variational form and that 
differential form is only formal and written for reader’s comfort. 

Effective dimensional equations: 

The model is given in terms of u®® : —>■ which is the vector of components of the 
displacement of the middle surface of the shell in the contravariant basis and p®® : —)• M 
which is the pressure in the 3D shell. Let us denote the bending moment (contact couple) 
due to the variation in pore pressure across the plate thickness by 

/3 2;/rv 

m = -C®(A®p(u®S))A® + / ysp^^dysA<^, (2.10) 

J-2 A + Ip J_£/2 


where p(-) is given by (| 2 . 6 I) and C® is the elasticity tensor, usually appearing in the classical 
shell theories, given by 

A 

_ e-r ( IhM I _L V//IH’. W. d A,, 

sym * 


C®E = 2p- 


tr (E)I + 2pE, E e 


\ + 2p 

Then the model in the differential formulation reads as follows: 


-{riafi + Y - Y = {'PYo. + {p/)a in wl, a = 1, 2, 

/3=1 V K=1 K=1 / 

Y ( - Y | = iV^^)3 + (Pl% in UJl, 

a, 0=1 \ K,= l ) 

1 2 

-{po^uf -h <9/3<®) - Y - Kpuf = 0 in UL, a, ^ = 1, 2, 

K,= l 


uf = 0,f = 1,2,3, 


du 


= 0 on duoL, for every t G (0,T), 


( 2 . 11 ) 
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( 2 . 12 ) 


f^G 


a 


A + 2//y dt 


dp' 


.eff 


— a 


2/i 


\ + 2pi 


/In®® h /l^Ti®® 

A-p(^)y3--^=0 


dt 


V d{yzf 


in (0,r) X (JL X (-£/2,£/2) 
kdp"^^ 


V dy3 

r)®® - 
P — PL,in 


= -Vl, on (0, T) XULX {{-i/2] U {£/2}), 


given at t = 0 . 


Here {'P^^)i,i = 1,2,3 are components of the contact force o r at in the covariant 
basis, Vl = Vl o X, = p 1 in ° Thus, the poroelastic flexural shell model in the dif¬ 
ferential formulation is given for unknowns {n, m, and by equations (| 2 . 10 l) . (| 2 . 11 l) 

and (|2.12l) . The components of n are the contact forces, linked to the constraint 7 (u®®) = 0, 
and being the Lagrange multipliers in the problem. The components of m are the contact 
couples. The first two equations in (I2.11h can be found in the differential equation of the 
Koiter shell model (see [3 Theorem 7.1-3]). The third equation is the restriction of approx¬ 
imate inextensibility of the shell. The first equation in (I2.12p is the evolution equation for 
the effective pressure with associated boundary and initial conditions in the remaining part 
of ()2.12l) . Note also that the same model holds for the shell clamped only on the portion of 
the boundary, i.e., the boundary condition in the fourth equation in (| 2 . 1 ip holds for a subset 
of dtOL with positive measure. 

In subsection 13.11 we present the dimensionless form of the problem. In subsection 13.21 
we recall existence and uniqueness result of the smooth solution for the starting problem. 
Subsection 13.31 is consecrated to the introduction of the problem in curvilinear coordinates 
and the rescaled problem, posed on the domain H = cj x (—1/2,1/2). In subsection 13.41 we 
formulate the main convergence results. In Section |4| we study the a priori estimates for the 
family of solutions. Then in Section [5| we study convergence of the solutions to the rescaled 
problem, as e ^ 0. In Appendix we give properties of the metric and curvature tensors. 


3 Problem setting in curvilinear coordinates and the main 
results 


3.1 Dimensionless equations 

We introduce the dimensionless unknowns and variable by setting 


/3 = /3G/i; P = ^-, Uu^ = u- T=^- A =-; 

L kpL pL 


ePU 


Pp^ = p; yL = y; x^L = 3 : 3 ; rL = r; XL = X; tT = t; 

1 j 


After dropping wiggles in the coordinates and in the time, the system ()2.7p - (l2.9p becomes 


- div cj^ = - A - A V div -I- a V = 0 in (0) T) x ( 3 . 1 ) 

= 2e(u^) -I- (A div — ap^)I in (0, T) x (3.2) 

d 

— {f3p^ + adw\V) — e‘^Ap'^ = 0 in (0,r)xH'^, (3.3) 


where W = (iif, hi, U 3 ) denotes the dimensionless displacement field and p^ the dimensionless 
pressure. We study a shell with thickness s = i/L and section cj = ul/L. It is described 
by 

O® = r({(xi,X2,X3) G w X (-e/ 2 ,e/ 2 )}) = r(H^), 
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(respectively SI) is the upper face (respectively the lower face) of the shell is the 

lateral boundary, = duj x (—e/2,e/2). 

We suppose that a given dimensionless traction force is applied on S^ U SI and impose 
the shell is clamped on 


= (2e(u^) — ap^I + A(div u^)/)i>' = s^V± on S^ U Si, 
= 0, on 


(3.4) 

(3.5) 


For the pressure p^, at the lateral boundary F^ we impose zero inflow/outflow flux: 

- \/p^ • 1/ = 0. (3.6) 

and at S^ U SI, we set 

-\/p^-u = ±V. (3.7) 

Finally, we need an initial condition for p^ at t = 0, 

p^(xi,X2,X3,0) = epin(xi,X2) m . (3.8) 

Let V(f2^) = {v £ : v|pe = 0}. Then the weak formulation corresponding to 

(I3T])-(I3SI) is given by 

Find £ H^{0,T,V{n^)), p^ £ such that it holds 




2 e(u^) : e(v) dx + \ 


/ div div Y dx — a div v dx 

J J 


k 


e^V+ • V ds + ■ v ds, for every v £ V’(f2^) and t £ (0,T), (3.9) 


/3 / dtp^qdx+ / a dtv q dx + / Vp^-Vqdx 

Jne Jqs Jqc 

= J Vqds — £^ [ Vq ds, for every qEH^{tl^) andt£(0,T), (3.10) 




p^|{t= 0 } =£Pm, in 


(3.11) 


Note that for two 3x3 matrices A and B the Frobenius scalar product is denoted by 
T : .B = tr {AB^). 


3.2 Existence and uniqueness for the £-problem 

In this subsection we recall the existence and uniqueness of a solution {u^, p^} eH^{0,T-,V{n^))x 
of the problem ()3.9l) - (l3.1ip . We follow [16] and get 

Proposition 1. Let us suppose 

Pin £ 77^(17^), V±e andV e H\0,T-,L‘^{uj)), = 0. (3.12) 

Then problem L3.9\) - I13.11\) has a unique solution {u^,p^} £ 77^(0, T; V(I7^)))x77^(0, T ; 77^(17^)). 


8 






3.3 Problem in Curvilinear Coordinates and the Scaled Problem 


Our goal is to find the limits of the solutions of problem (I3.9I) - (I3.11I) when e tends to zero. 
It is known from similar considerations made for classical shells that asymptotic behavior 
of the longitudinal and transverse displacements of the elastic body is different. The same 
effect is expected in the present setting. Therefore we need to consider asymptotic behavior 
of the local components of the displacement u*". It can be done in many ways, but in 
order to preserve some important properties of bilinear forms, such as positive definiteness 
and symmetry, we rewrite the equations in curvilinear coordinates defined by r. Then we 
formulate equivalent problems posed on the domain independent of £. 

The covariant basis of the shell which is the three-dimensional manifold parameterized 
by r, is defined by 

gf = diV : M^, i = 1, 2,3. 

Vectors {gf,g 2 ,g 3 } are given by 


51 = My) + xsdyi^siv), 

52 = My) + xsdy^Mv), 
g| = My)- 

Vectors {g^’*',g^’^,g^’^} satisfying 

-eI = 6ij on Ti\ i,j = 1,2,3, 


where 5ij is the Kronecker symbol, form the contravariant basis on P . The contravariant 
metric tensor = {g^^’^), the covariant metric tensor = Mj) and the Christoffel 

symbols T*’^ of the shell P are defined by 


g^3,e ^ 


g 




ah = gf 


g 


• djgl on P , i,j,k = 1,2, 3. 


We set ^ 

r*’" = (rJ’D, >=!,...,3 and 7 ^(v) = i(Vv + Vv^)-^u,r’^ (3.13) 

i=l 

Let = det G^. Until now we were using the canonical basis { 61 , 62 , 63 }, for M^. Now the 
displacement is rewritten in the contravariant basis. 


33 3 

u^or{yi,y 2 ,X 3 ) = o r{yi,y 2 ,X 3 )ei = ^ <( 2 / 1 , y2, a;3)g*’^(l/i, 2/2, 3 : 3 ), vor = 

i=l i=l i=l 

while for scalar fields we just change the coordinates 


p^or=p^, qor = q, V o r = V, Pin o r = pm, 

on . The contact forces are rewritten in the covariant basis of the shell 

3 

V±oY = '^{V±)igl on 
2=1 


New vector functions are defined by 


W = uf6j, V = ViGi, V± = {V±)iei. 
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Note that uf are not components of the physical displacement. They are just intermediate 
functions which will be used to reconstruct u^. The corresponding function space to V(n^) 
is the space 

= {v e : v|re = 0}. 

Let = (Vr)”"^ and let 

CE = A(trE)I + 2E, for all E G (3.14) 

Then the problem (I3.9p - ()3.1ip can be written by 

[ C (Q^ 7 ^(u^)(Q^)^) : (Q^7£(v)(Q^)^) V¥dy -a [ /tr(Q^ 7 ^(v)(Q^)^) 

= V+■v^/^ds + s^ / V-^w^/^ds, V G V(n^), a.e. tG[0,T], 

dsi 

Jne a^tr(Q^7£(u'')(Q^)^)yv^dy 

+ • Q'^Vqy/^dy = [ Vqyf^ds - [ Vqyf^ds, 

dne Jsi 

q G a.e. t G [0,r], 

= Spin, for t = 0. 

(3.15) 

This is the problem in curvilinear coordinates. 

Problems for all u^,p^ and are posed on e-dependent domains. In the sequel we 

follow the idea from Ciarlet, Destuynder [8] and rewrite (j3.15p on the canonical domain 
independent of e. As a consequence, the coefficients of the resulting weak formulation will 
depend on e explicitly. 

Let n = cu X (—1/2,1/2) and let : fl —)• be defined by 

K^{z) = (z^,z^,ez^), z G Q, e G (0,eo). 


By = u! X {±1/2} we denote the upper and lower face of fl. Let T = doj x (—1/2,1/2). 

To the functions u^, p^, , gf, g^’*, Q'^, T*’/, i,j,k = 1,2,3 dehned on fl we associate 
the functions u(e), p{e), y(e), gi(e), g*(e), Q(£), hj,k = 1,2,3 defined on Q by 

composition with R^. Let us also define 

V(n) = {v = (ui,U 2 ,U 3 ) G R^(0;M^) : v|r = 0}. 

Then the problem (13.151) can be written as 


C (Q(e)7^(u(e))Q(e)^) : (Q(e)7''(v)Q(e)^) y/^dz 


— ea 


[ p(e)tr(Q(e)7^(v)Q(e)^)\/®d^ 

Jn 


= • VY^y(eyds ± • vy^ g{s)ds, v G V’(n), a.e. t G [0, T], 

Jt,+ Jt,- 

e / f o;^tT{Q{e)Y{u{£))Q{£f)q^/^dz 

J J 

±6^/ Q(e)V^p(e) • Q,{£)V^q^g{£)dz 
Jn 

= / Vq\fg{^ds — £^ [ Vqyjg{£)ds, q G a.e. tG[0,T], 

4s_ is+ 

p{£) = £Pin, for t = 0. 


(3.16) 
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Here 


7^(v) = \iz(y)+iyiy) - 

i=l 


(3.17) 



0 

0 

^dsvi ' 


divi 

^{d2Vi +diV2) 

\dlv3 ‘ 

7 z(v) = 

0 

0 


, 7 j;(v) = 

^{d2Vi +diV2) 

d2V2 

\d2v3 


_ I53U1 I93U2 

dsvs 


\d2v3 

0 


+Vj^g, = [ 0 0 dsq], Vyq=[diq d2q O]. 

We assume for simplicity that pi^ = 0 and proceed with asymptotic analysis, 
by rescaling the pressure 




We start 


The equations are now 


e [ C (Q(e) 7 ^(u(e))Q(e)^) : (Q(e) 7 ^(v)Q(e)'^) ^/^dz 
Jn 

-e^a [ 7r(e)tr(Q(e)7''(v)Q(e)^)v^^(iz 
Jn 

= / V± ■ 'v^Jg{e)ds, v G V(H), a.e. t G [0,T], 

Jt,± 

/ P^^^qV^dz + e^ f a^tr{Q{e)Y{u{e))Q{s)^)q^/^dz 

t/ i/ 

+ e^ [ Q(e)V^ 7 r(e) • Q(e)V^gv^^d 2 ; = 76 ^ / Vq^fg{e)ds, 
J Q, J S-i- 

q G a.e. t G [0,T], 

p(e) = 0, for t = 0. 


(3.18) 


Here and in the sequel we use the notation 

T / Vq^/g{£)ds = / Vq^/g{e)ds - / Vq^/g{£)ds, 

Jt.± Jt,- Jt.+ 

[ 'P± ■ M^/^ds = f V+- ^r^/^ds + f V-- ^^/^ds. 

JYi± JYaJ^ JY,— 


Remark 2. Existence and uniqueness of a smooth solution to problem (I3.18|) follows from 
Proposition [1] and the smoothness of the curvilinear coordinates transformation. 


3.4 Convergence results 

In the remainder of the paper we make the following assumptions 

Assumption 3. For simplicity, we assume thatpin = 0, thatV G H^{0,T] L'^{lo)), = 

0 and thatV± G 77^(0,T;L^(a;;M^)), withV±\^^^Qy = 0. 

We recall that the differential operators 7 and p are given by (12.51) and ()2.6p . respectively. 
Let 

Vf(w) = {v G H^iuj) X H^{uj) X H^{uj) : 7 (v) = 0, v|a^ = 0, -^\duj = 0}. (3.19) 

We will suppose the classical hypothesis that leads to the ’’flexural shell” models: 

Vf{uj) / {0}. (3.20) 
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Let us formulate the boundary value problem in O = cu x (—1/2,1/2) for the effective dis¬ 
placement and the effective pressure: 

Find {u, vr*^} G (^([O,T]; Vir(w) x L^(ll)), G L^((0,T) x fl) satisfying the system 


1 

12 


C(A‘^p(u)) : p(v)A'^-v/odzid ^;2 + 


2q; 


1/2 


= / (V+ -I- V-) • vy/adzidz 2 , v G Vir(a;). 
J (jJ 


zsTT^dz^A'^ : p{'v)y/adzidz2 
\ + 2JuJ-i/2 ^321) 


s. L + 


a 


\ + 2 


IT 


) /-, f 2a du f dn^ dq /- 

qy/adz- / ~-A ■. p{ — )z^qy/adz + / -^—-^^adz 

Jfi A + 2 dt Jq dz3 dz3 


= T 


f Vqy/a dzidz 2 in T>'{0,T), q € H^{—1/2,1/2;L‘^{lo)), 

Jt,± 


TT^ = 0 at t = 0, 


(3.22) 


(3.23) 


where 7 (-) and p(-) are given by (|2.5p and (|2.6I) . respectively, and 


CE = 2^^ tr (E)I + 2E, Eg (3.24) 

X + 2 ^ 

Proposition 4. Under Assumption\^ problem \3.21\) ~ \3.23\) has a unique solution {u, vr*^} G 
C{[0,T]-,Vf{<^) X L?{yt)), 5^37r° G L^((0,T) x Ul) Furthermore, dtir^ G L^((0,T) x 11) and 
9tuGL2(0,r;Vir(w))). 


Proof. First we prove that {u, vr*^} G C([0, T]; Vf(^^) x L^(11)) and G L^((0, T) x 11) 

imply a higher regularity in time: 

Let us take q = Z 3 q{zi,Z 2 ), q G C°°iuj) as a test function in (|3.22l) . It yields 

^ dz 3 - : p(u) G H\0,T; L\u;)). (3.25) 

V A + 2/ J-i/2 6 A + 2 

After inserting (13.251) into (I3.2ip . it takes the form 


f C(A'^p(u)) : p{^/)A'^^/adzldz 2 ^ / A'^ : p(u)A^ : p(y)y/adzidz 2 

Ju) ^ Juj 

= / {V++7^^) ■ \'y/adzidz 2 — (^ / F{t)A‘^ : p{v)^/adzldz 2 , v G Vf(w), 

J LJ J (jJ 


(3.26) 


with = a^/(A + 2)/(/3(A-|-2) and c| = 2a/{l3{\ + 2) + a^). Taking the time derivative 
and using the time regularity of F and V±, yields dtu G L^(0, T ; Vir(a;)). For such u classical 
regularity theory for the second order linear parabolic equations applied at (I3.26P implies 
dtTT^ G L2((0,r) X H). 

clu 

The existence and the uniqueness are based on the energy estimate. If we choose v = 

ot 

as a test function for equation (13.211) and vr^ as a test function in (I3.22p and sum up the 
equations to obtain the equality 



/ dt{V+ + P_) • u^/adzldz2 =F / Vir^y/a dzidz2- 

J uj J S-1- 


(3.27) 
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Equality (j3.27l) implies uniqueness of solutions to problem (I3.21I) - (I3.‘23|) . Concerning exis¬ 
tence, equality (j3.27p allows to obtain the uniform bounds for p(u) in for 

7 r° in L°°{0,T; L‘^{Q)) and for in L‘^{0,T; L‘^{Q)). Using [51 Teorem 4.3-4.] and the 

classical weak compactness reasoning, we conclude the existence of at least one solution. □ 

Remark 5. Let (3 = (5 + Using separation of variables, we obtain the formulas 


4 f—1)-^ 


* 7r^(2j — l)^(t — r), d . . 

exp{-=- — zi,Z 2 ) 


2 a 

\ + 2 

0 . U 8 

Z^TT dZ3 = + 


13 ' dr 

'Zot \ 

+ r--A'" : p(u(r,zi, 2 : 2 ))) dr] sin((2j - l)7rz3) 

\ —I— A ' 


(3.28) 


/- 1/2 


12 


TT 


I +°o .. 


7 r^( 2 j — l) 2 (t — r) d 

— }^{PV{T,ZUZ2) 




2 a 


+ r --A'" ; p(u(r, 21 , 2 : 2 ))) dr. 

A + 2 


(3.29) 


After plugging formula \3.21A) into equation U3.21\) . we observe memory effects in the flexion 
equation. 

The main result of the paper is the following theorem. 

Theorem 6. Let us suppose Assumption\M Tet {u(e), 7r(e)} & H^{0,T;V{Ll))xH^{0,T-, H^{Ll)) 
be the unique solution of 113.18\) and let {u, vri'} be the unique solution for \3.21\) - \3.23\) . Then 
we obtain 

u(e) —)• u strongly in (^([O, T]; M^)), 

- 7 ^(u(e)) —)• strongly in (^([O,Tj;L2(n;]R^^^)), 


7r(e) —7r° 
97r(e) diT^^ 


where 


dz^ 


7° = 


dz?, 


strongly in C{[0,T];L‘^{LI)), 
strongly in L‘^{0,T; L‘^{Q)), 


-23P(u) 
0 0 


0 

0 

^ 7 r° 23 ^rA_A'^ : p(u) 


\ + 2 


\ + 2 


(3.30) 


As a consequence of the convergence of the term - 7 ^(u(e)), we obtain the convergence of 
the scaled stress tensor. 


Corollary 7. For the stress tensor a(e) = C(Q(e)j^(u(e))Q(e)'^) — ap(e)I one has 
1 


( 7 (e) —)• (7 = C(Q 7 °Q^) — aF^l strongly in (^([O, T]] L"‘. 


3x3i 


(3.31) 


The limit stress in the local contravariant basis Q = (a^ a^) is given by 


Qflaq, = 


2 a 

\ + 2 


- 23 


2A 

\ + 2 


(A^ : p(u))I + 2A"p(u) A'^ 0 


0 


0 
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4 A priori estimates 


Fundamental for a priori estimates for thin shell-like bodies is the following three-dimensional 
inequality of Korn’s type for a family of linearly elastic shells. 

Theorem 8 ([71 Theorem 5.3-1], (TOl Theorem 4.1]). Assume that X e Then 

there exist constants eo > 0, C > 0 such that for all e € (0,eo) one has 

(j 

I|v11h 1(Q;R3) < —ll7^(v)lli2(Q.R3x3), Vv G V(n). 

Remark 9. Only a portion of the boundary with positive surface has to be clamped for the 
statement of the theorem to hold. 


Now we state the asymptotic properties of the coefficients in the equation (I3.18p . Direct 
calculation shows that there are constants rUg, Mg, independent of e G (0,eo)j such that for 
all z G D, 

mg < < Mg. (4.1) 

The functions gi*(e),gij(e), 5 ®'^ (e), 5 (e), Q(e),r*^(e) are in (7(0) by assumptions. Moreover, 
there is a constant (7 > 0 such that for all e G (0,eo)) 


\\g\e) - a*lloo + \\giie) - ctilloo < Ce, 
\\-^^/^\\oo + ||\/5(^- \/«l|oo < Cs, 


(4.2) 


where 


llQ(e)-Q-ez3Qilloo<C£2, = ( ELi ELi O), 

Iji (rj,(e) - r;.,(0)) - ^3 WH- < Ce, 

is the norm in (7(0). For proof see |9]. Additionally, in [71 Theorem 3.3-1] the 


asymptotic behaviour of the Christoffel symbols in L°°-norm is given by 


T^{e) = 


Ff^ - szsbfh 

T^^-EZ3bf\2 
2 

-bf - ezs'^blbf -b^ - ezs'^blb 

T = 1 

where k = 1, 2 and 

6 ii- 


V 


Tf^-ezab^h -bf - ezs'^bf 

r=l 

2 

T ^2 - b -^2 - ^^3 bib 


T=1 


T = 1 


+ 0(E), (4.3) 


r'(e) = 


£Z3 Y ^12 - ^^3 Y 0 


K =1 

2 


K =1 

2 


b2l - £Z3 Y ^2bKl b22 - £Z3 Y ^ 2^«^2 0 

K=1 K=1 


(4.4) 


0 0 0 
In the following two lemmas we derive the a priori estimates in a classical way. 
Lemma 10 . There is C > 0 and Eq > 0 such that for all e G (0,eo) one has 


,1 


7^(u(e))llL°°(0,T;L2(Q;R3x3)), lk(e) ]] 2,00 (o,T;L 2(0;R)) i IIl2(0,T;L2(Q;R3)) < C. 
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{r') 

Proof. We set v = —and g = 7r(e) in (I3.18P and sum up the equations. After noticing 
that the pressure term from the first equation cancels with the compression term from the 
second equation we obtain 

(Q(e)7^(u(e))Q(e)^) : (Q(e) 7 "(u(e))Q(e)^) yf^dz 

+ X In (4.5) 

= 6^/ V± • y/g{e)ds p / Pvr(e)v^c/(e)ds. 

J J S-|- 

Dividing the equation by and using the product rule for derivatives with respect to time 
on the right hand side we obtain 

(XX*^ (Q(e)7^(u(e))Q(e)^) : (Q(e) 7 "(u(e))Q(e)^) ^/^dz + P ^/^dz'^ 

+ e‘^ [ Q(s)V^'k{s) ■ Cl{s)V^7r{s)-s/g{s)dz 

Jn 

= 4 /" • u(e)v^5(e)(is - f ■ u{£)^/^ds^ f V^T{£)^/^ds. 

Now we use the Newton-Leibniz formula for the right hand side terms and the notation 

V = {V++ V-)Z3 + ~ , V = 2Pz3 


to obtain 

(XX^ (Q(e) 7 ^(u(e))Q(e)'^) : (Q(e)7^(u(e))Q(e)^) y/^dz +/3 J^ 7 r{e)^ y/^dz'^ 

+ £^ [ Q(e)V^ 7 r(e) • Q{e)V^ 7 r{e)y'g{£)dz 
Jn 

= 5 X ■"(")-1 ^ ■ “(")v^) 

Next we integrate this equality over time 

^4 [ (Q(e) 7 ^(u(e))Q(e)'^) : {Q{£)'Y%u{£))Q{£f)^/^dz + ^|3 f irief^'^dz 

Jn ^ Jn 

+ / / Q(e)V^ 7 r(e) • Q(e)V^ 7 r(e)-y/ 5 (e)(izdr 

Jo Jn 

= ?f 4 / (Q(^) 7 ^(u(e)|t=o)Q(e)'^) : (Q(e) 7 ^(u(e)|t=o)Q(e)'^) +/3 [ T^{£f\t=o\/^dz 

+ [ 4-(^ • u(e)v^^)(iz - [ -^{V\t=o ■ u{e)\t=o^/^)dz 
Jn <2^3 20 <22:3 

-ULirX^' - i‘ X ^ d.ir. 

(4.6) 
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Since we have enough regularity for u(e) we consider (I3.18P for t = 0. Then u(e)|i=o satisfies: 
for all V G V(r2) 


4 / C (Q(e)7^(u(e)|t=o)Q(e)^) : (Q(e) 7 ^(v)Q(e)'^) y/^dz 
e Jn 

--a f 7r(e)|t=otr (Q(e) 7 ^(v)Q(e)^)v^^(iz = [ V±\t=o ■ v\/^ds. 
£ Jn Jt,± 


Since the initial condition is 7r(e)|t=o = 0 this equation is a classical 3D equation of shell-like 
body in curvilinear coordinates rescaled on the canonical domain. Next, V±\t=o = 0 and the 
classical theory (see Ciarlet [7]) yields u(e)|t=o = 0. Using Korn’s inequality, positivity of C 
and uniform positivity of Q(e)^Q(e) and g{e) in (I4.6p yields the estimate 

(Q(e)7^(u(e))Q(e)^) : {Q{e)Y{u{e))Q{ef)^/^dz + ^/3 j ^{ef^/^dz 

+ [ f Q(£)V^7r(e) • Q{£)V^7r{e)y^g{e)dzdT < C. 

io 4c 

Since C is positive definite and since g{£) is uniformly positive dehnite (see [71 Theorem 3.3-1]) 
we obtain the following uniform bounds 

“IIQ(£) 7 ^(u(e))Q(e)^||Loo(o 7 ;L 2 (f^.]R 3 x 3 )), ||vr(e)||x,oo(o7;i,2(Q)), ||eQ(e)V^ 7 r(e) 112,2(0772(^^.^3)). 

Since Q(e)^Q(e) is uniformly positive definite these estimates imply uniform bounds for 

-|| 7 ^(u(e))Q(e)^|| 2 ,oo(o 772 (f 2 .R 3 x 3 )) and lkV^ 7 r(e) 112 , 2 ( 0772 (^.^ 3 )). 

Applying the uniform bounds for Q(e)^Q(e) once again implies the statement of the lemma. 

□ 


We now first take the time derivative of the hrst equation in p3.18p and then insert 

V = —as a test function. Then we take q = as a test function in the second 

dt 

equation in (|3.18p and sum up the equations. The following equality holds 

+ I 3 J / +Q(e)V^7r(e) • Q(e)V=7r(e)x/^dz (4.7) 

Similarly as in Lemma [in] from this equality we obtain 

Lemma 11. There is C > 0 and eo > 0 such that for all e G (0,eo) one has 

ll-V7£_£|l 

ll “7 f—^^j||L2(o,r72(C;R3x3)), ||—^^||2,2 (o 772(0;R)), ||eV vr ||Loo(o772(f2;R3)) S L. 

As a consequence of the scaled Korn’s inequality from Theorem | 8 | we obtain 
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Corollary 12. Let us suppose Assumption\^ and let {u(e),7r(e)} be the solution for problem 
113.18\) . Then there is C > 0 and eo > 0 such that for all e G (0,eo) one has 


||-7^(u(e))||Hl(0,T;L2(Q;R3x3)), ||u(e)||Hl(0,T;/fl(O;R3)); (0,T;L2(f7;K)), 

II^^IIl°°(0.T;L2(O;1R)) S g. 

Furthermore, there are u G H^{0,T-, G T;L^(n;M)) and 7 ° G 

i/^(0, T; L^(Q; such that on a subsequence one has 

u(e) ^ u weakly in T;M^)), 

7r(e) ^ TT^ weakly in 17^(0, T; L^(n; M)), 

^ ^ ~Q — weakly in LF‘{t),T-,LF‘{Ll]'E3)) and weak * in L°°(0, T; M)), 
^ 7 ^(u(e)) ^ 7 ° weakly in 

Proof. Straightforward. □ 


Since i 7 ^(u(e)) depends on u(e) one expects that the limits u and 7 *^ are related. The 
following theorem gives the precise relationship. It is fundamental for obtaining the limit 
model in the classical flexural shell derivation as well as in the present derivation, see |10j . 
The tensor 7 is the linearized change of metric tensor and p is linearized change of curvature 
tensor. They usually appear in shell theories as strain tensors. 

Theorem 13 ((Tj Theorem 5.2- 2 ], m Lemma 3.3]). For any v G V(n), let^^{v) G 
and let the tensors ^{w), p{w) belong to , H~^ , respectively. Let the fam¬ 

ily {w(e)}£>o C V(n) satisfies 

w(e) ^ w weakly in 

- 7 ^(w(e)) ^ 7 weakly in 

as £ 0. Then the limit function w is independent of transverse variable Z 3 , belongs to 

(w) X X satisfies the clamping boundary conditions 


wlat,; = 0, 


and the following conditions 


dw 3 

du 


Iduj — 0 


dja3 


7 (w) = 0, p(w) G L (0;M ^ ) and —= -pafii^) 

0 Z 3 


If in addition there is x ^ such that as £ ^ 0 


p(w(e)) —X strongly in H 


then 

w{£) ^ w strongly in and p(w) = x £ 
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Remark 14. The estimates from Lemma [TU] and Lemma [TT] yield uniform boundedness of 
u(e) inC'°’^/^([0, r],V(fl)), i 7 '^(u(e)) in T]; and 7 r(e) in T]; L^(n)). 

Hence by Corollary [T^] and Aubin-Lions lemma (see [21]1. there is a subsequence such that 
the {u(e)} converges to u also in (^([O, T]; L^(n; M^)). 

Let if G L^(H). Then for every <5 > 0, there exists cps £ such that \\(p — p> 5 \\L 2 (n) < 

6 . Next 

f ^ 

o<t<T Jn Oxi 

< sup I /" ^(u(e)(t) - u(t))(v 9 - 995 ) dx| + sup | [ ^^{u{e){t) - u{t)) dx\ (4.9) 

0 <t<T Jq UXi 0 <t<T Jn UXi 

< C'5||u(e) - u||c([o,r];^^bf^;K^)) 
for £ < £o{ 6 ). Therefore 

r Q _ 

lim sup I / TT—(u(e)(t) — u(t ))(/7 dx\ < C 6 , 

^^^o<t<T Jn 

which yields 

u(e)(f) ^ u(f) weakly in R^(H;M^) for every t G [0,T]. (4.10) 

Argument for the sequences {- 7 ^(u(e))} and {vr(e)} is analogous and we get 
7r(e)(t) ^ weakly in L^(H), 

i 7 =(u(e))(t) ^ 7 °(t) weakly in 1 . 2 ( 11 ;M^xs^ 

for every t G [0,T]. 

Thus we may apply Theorem [T3j with w(e) = u(e)(t), for each t G [0,T] and conclude 
that the limit points of {u(e)(t)} belong to Vf(w)- 
Moreover we conclude that 

Ta/3 Tck/3 •2^3Pq.^(u), 

where 7 ^,^ do not depend on Z 3 . We denote 7 = Yia 0 \a,ii=i, 2 - 

5 Derivation of the limit model 

In this section we derive a two-dimensional model. We obtain it in five steps. In the first 
two we take the limit in (I3.18P for special choices of test function. In this way we addition¬ 
ally specify the limits u, 7 r^, 7 ^ and the equations they satisfy. The part 7 of 7 ^, which is 
independent of 2 : 3 , is identified in Step 3 by techniques usually applied in the proof of strong 
convergence of strain tensors in the classic shell models derivations. In Step 4 we prove the 
strong convergence of displacements, while in Step 5 we prove the strong convergence of stress 
tensors. 

Step 1 (Identification of 7 ^ 3 ). Now we are in a position to take the limit as e ^ 0 in 
(j3.18l) with the first equation divided by e: 

j^C ^Q(e)^ 7 ^(u(e))Q(e)'^^ : (Q(e)e 7 ^(v)Q(e)'^) y/^dz 

-a [ 7r(e)tr (Q(e)e 7 =(v)Q(e)^)V^ 5 (e)d 2 ; = ^2 / V±■v^/^ds, 

Jq J Sj- 

V G V(H),t G [0,T]. 
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In the limit we obtain 


C 


(Q( 0 ) 7 °Q( 0 )'^) : (Q(0)72(v)Q(0)'^) ^/^dz -a [ 7 r°tr (Q( 0 ) 7 ^(v)Q( 0 )'^) = 0, 

Jn 


v€V{n),t€ [0,T], 

which, nsing Q(0) = Q and g{0) = a, yields 

[ [C (Q 7 °Q^) — avr^l) : (Q 7 j,(v)Q'^) y/adz = 0, v e V(n), a.e. t € [0,r]. 
Jn 

From the definition of 7 ^ and the function space V(n) we obtain 

(Q^ (C (Q7°Q^) - aTfOl) Q)^3 = 0, z = 1,2,3. 

This implies 

'a tr (Q 7 °Q^) - a 7 r°) Q^Q + 2 Q^Q 70 Q^q) =0, z = 1, 2, 3. 


i3 


Since 


Q^Q = 


A‘= 0 
0 1 


(5.1) 


we obtain expressions for the third column of 7 *^ in terms of the remaining elements 

(Q^Q 7 °)i 3 = (Q^Q7°)23 = Atr (Q'"Q 7 °) - an^ + 2^1^ = 0 . (5.2) 

The first two equations imply that 


7i3 

723 J 


= 0 


and since is positive definite we obtain that 75^3 = 733 = 723 = 732 = 0. Prom the third 
equation in ()5.2jl we obtain 


AA7 


711 7i2 

712 722 


— avr^ + (A + 2)733 — 0- 


Thus we have obtained the following result. 

Lemma 15. 

7?3 = 731 = 723 = 732 = 0> 
A 


..,0 _ “ „o 

733 — ^^ “ 


-A7 


7?i 7?2 
7 i2 722 


A+2 A+2 

From this lemma and Theorem n we have that 7 ^ is of the following form 


7 O = 


7 - Z3p{^) 
0 0 


0 

0 

jf^vr0-^A-(7-.3P(u)) 


Step 2 (Taking the second limit). Let v e Vf(w), where Vf(w) is given in (j3.19p . and 
let be given by 

vl{z) = -{diV 3 + 2vib\ + 2 v 2 bl)z 3 , 
vliz) = -{d2V3 + 2vib\ + 2v2bl)z3, 
vl{z) = 0. 
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Then 


7^(v) = 7z(v^) + 7j/(v) - ^nir*(0) = 0 


2=1 


and v(e) = v + e V(r2). 

A simple calculation shows that ^ 7 '"(v(e)) = ©(v) +eG, where 

®(v) = 7j/(v^) - r*(0) - (0)) 

and G is bounded in L°°{U]'R^^^). For the test function v(e) equation ()3.18p now reads 
J^C ^Q(e)^7"(u(e))Q(e)'^^ : (^Q(e)^7^(v(e))Q(e)^^ \/W)dz 


- a ^ 7 r(e) tr ^Q(e)^7=(v(e))Q(e)'^^ \/W)d. 


z = 




'P± ■ v(e)\/ff(e)(is, 


Jj?^<,sMZ)dz 


d 


m 


dt 


1 


Q(e)-7^(u(e))Q(e) 



£]dz 


+ e^ / Q(e)V^7r(e) • Q{£)V^qy/g{£)dz = =F / Vq^/g{£)ds, 

Jn Jt.± 

q S a.e. t G [0,r]. 

In the limit when e ^ 0 we obtain 

f C (Q7 °Q^) : (Q0(v)Q'^) y/adz — a / vr^tr (Q0(v)Q^) 

J Q J Q 

= / {V++V-) ■ \-y/ads, V G d/g(ca;M^), a.e. tG[0,r], 

J u) 

J^f3^^qy/ddz + j a-^tr[Q'y^Q^)qy/ddz + j ■-^Qesy/ddz 

= T [ VqVdds, q e H\n),t e [0,T]. 

Jt.+ 


(5.3) 


Note that Qeg • Qeg = 1. 

According to Lemma [20] (in the Appendix) one has 


©(v) = -Z3 


P(v) 




'^b^{dKV3 + '^Vrbl) 

K,= l T=1 

2 2 


n=l 


T = 1 


K,= l 


T=1 


K=1 


T = 1 


Thus, using (|5.1D we obtain 

tr(Q©(v)Q^) =tr(Q^Q©(v)) = -Z3A7p(v). 


(5.4) 
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Next, using Lemma [TS] and Remark 1141 we compute 
tr(Q7°Q'^)=tr(Q^Q7°) = A'^: 

= A": 

2 


7ii 7 i2 
L 7i2 722 


+ 733 


711 7i2 

712 722 


+ 


—ttO - ^A7 

A+2 A+2 


711 7i2 

712 722 


A + 2 
2 

A + 2' 


A7 

A7 


711 7i2 

712 722 

711 7i2 

712 722 


+ 


a 


-TT 


A + 2 
2 

^^A + 2 


A7p(u) + 


a 


A + 2 


-TT . 


Further, using (15.ip and Lemma [15] we obtain 

Q^Q7°Q^Q = 


A''(7 - ^3P(u))A'= 0 

0 


A^(7-^3P(u))A'^ 

0 


733 


A+2 A+2 


A+(7-^3P(u)) 


The main elastic term is now 


f C (Q7°Q^) : (Q0(v)Q^) y/adz 

Jn 

= / Atr (Q7°Q'^) tr (Q©(v)Q'^) + 2Q7°Q^:Q©(v)Q'^v^(iz 

Jn 

{-Z 3 A'^:p{v)) y/adz 


= f A f ~-^A'=:7 - Z3^^A7p(u) + “ 

Jn 


n V A + 2 


A + 2 


A + 2 


+ [ 2Q^Q7°Q^Q:©(v)V^dz 
Jn 

= f ( 23 )"t^A=:p(u)A‘:p(v) + -» 

Jn 


A + 2 


A + 2 


7r° (-Z 3 A'^:p(v)) y/adz 


+ f 2{z3)‘^A^p{u)A'^:p{v)y/adz. 

Jn 

Using the tensor C, defined by (|3.24p . the elastic term can now be written by 

[ C (Q7°Q^) : (Q©(v)Q^) V^dz 

Jn 


f C{A^p{u)) : p{v)A‘'y/ddzidz 2 + [ 
Ju) Jn 


Xa 

f2 A + 2 


7r° (- 23 A'^:p(v)) y/adz. 


The first equation from (15.31) now becomes: for all v £ V+Po;) one has 

-7r° (— 23 A'^ : p(v)) y/adz 


^ [ C{A''p{u)) : p(y)A‘'y/ddzidz 2 + ^ _o 


/f2 A + 2 


— a 


i; / 7r° {—z^A^ : p(v)) y/adz = / (R+ +V-) ■ \-y/adzidz 2 . 

kJ kJ (jJ 

This implies 

^ j C{A''p{u)):p{v)A''y/ddzidz- 


Z2 + 


2q; 


+/2 


= / (R+ + R_) • \^y/adzidz 2 , v £ Vf{u}). 
J LJ 


I . , zsTT^dzsj A^:p{-v)y/adzidz2 

X + 2 Juj \ J-l/2 J 


(5.5) 


(5.6) 


(5.7) 


(5.8) 
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The matrix 7 does not appear in (15.81) and is not important in the classical shell theory. 
However this is not the case here since the term 7 appears in the second equation in (|5.3I) . 
It will turn out to be 0 in the proof of the strong convergence, see Step 3. The equation (15.8p 
appears in the classical flexural shell model without the pressure term. 

The second equation in (|5.3I) can be now written as 


/ 

JO 


/3^—q^/adz+ / a 


dt 


d 


dt\X + 2 


■ 'Y - 


A + 2 


: p(u) + 


a 


A + 2 


7 r° ) qy/adz 


+ [ ^^-^y/adz = ^ f Vqy/ads, qeH^{fl). 
Jq dz3 dz3 


Then 


/ 

Jo 


+ \ ^^0 


„ a“ \ O'K'^ . f d 

P + - -zrr(lyo.dz + / a— . ~ 

x + 2) dtyx +2 


- 2:3 - ^ ^ A^:p(u) ) qy/adz 


A + 2 


+ [ t;— -^^y/adz = + / Vqy/ads, q € 
Jq dz3 dz3 


diT^ dq 


(5.9) 


A flexural poroelastic shell model will follow from (15.81) . (15.91) once 7 is determined. 
Remark 16. Let us set 


= Tr^ + 


2a 


/3(A + 2 ) + o? 


A" :7 


Then the couple {u, P^} satisfies the system 


— / C(A'^p(u)) : p{v)A'^y/adzidz 2 + 


2 a 


12 


r.1/2 


Z3P^dz3 ) A'^ : p(y)y/adzidz2 


+ 


A + 2 Ju \ J—1/2 

= / (P+ + P-) ■ yy/adzidz 2 , v G Vf(w). 
J U) 

2 \ dP^ f 2 _ 5 


a 


A+ 2 / dt 
dP^ dq 


f 

-qy/adz — / 
Jo 


azs^-A'^ : —p{u)qy/ddz 

o A + 2 dt 


„ „ y/adz = + / Vqy/ads, q G 

Q dz3 dz3 


(5.10) 


(5.11) 


Analogously to Section \3.4\ we prove that system 115.10\) - I15.11\) has a unique solution. It 
yields convergence of the whole sequence {u(e)}. Unfortunately, it is still not enough to have 
conclusions for {7r(e)}. 

Step 3 (Identification of 7 the strong convergence of the strain tensor and the 
pressure). We start with 

(^7^(u(e))(t) -7°(t)) Q(e)^) 

: (^Q(e) Q7^(u(e))(t) - 7°(i)^ Q(e)'^^ ^(7r(e)(t) - TT°{t))‘^y/^dz 

+ e^ f f (V"7r(e) - V"7rO)Q(e)^ • (V"7r(e) - V^ 7 r^)Q{e f y^dz. 

Jo Jo 

We will show that A(e) —>• A as e tends to zero. Since A(e) > 0 the A > 0 as well. After 
some calculation we will show that actually A = 0 and 7 = 0. This will imply the strong 
convergence in (|4.8p . 
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Since we have only weak convergences in (14.81) we first remove quadratic terms in A(e) 
using ()4.5p divided by Integration of ()4.5p over time, using 7r(e)|t=o = 0 and u(e)|t=o = 0, 
implies 

(^Q(e)i7^(u(e))Q(e)^^ : (^Q(e)i7^(u(e))Q(e)^^ \fg{e)dz 

+ / TT{£f ^/^dz + [ [ Q{£)V^7r{£) ■ Q{e)V^^^{e)^/^dzdT 

z Jn Jo Jn 



V± 


0 Jt,± 


■ ^^^V^dsdT 


0 Jn 

=F / / VTr{£) y/giFjdsdr. 
J 0 


Inserting this into the dehnition of A(e) we obtain 


K{e){t) = 


0 JT,± 


'P± • V g{e)dsdT 7 [ [ V7r{£)y/g{e)dsdT 

1 

I-' 

e 


0 JT:± 


- J^c (^Q(e)- 7 ^(u(e))(t)Q(e)^ ) : (Q(e) 7 °(t)Q(e)'^) yf^dz 

-/3 [ 'K{£)'n^\/g{£)dz 
Jn 


- 2e" 


I 



0 Jn 


Q(e)V^7r(e)(t) • Q{£)V^TT^{t)y/g{£)dzdT 


1 


+ 2 / ^ (Q(e)7°Q(e)^) : (Q(e)7°Q(e)^) 'JWidz + / i'^^fVWidz 

u t/ 


-he" 


Q(e)V^7r° • Q(e)V^7r°y^ g{£)dzdT. 


0 Jn 


Now we take the limit as e tends to zero and obtain that A(e) —>• A > 0, where 


+ V-) ■ %-^/adsdT 7 
at 


V y/adsdr 


0 JT.± 


t; [ C (Q7°Q^) : (Q7°Q^) Vadz - ;^/3 / {Tr^fy/adz - [ [ 
z Jn z Jn Jo Jn 


SvrO 

9zs 


(5.12) 


y/adzdr. 


d\i 


We now insert as a test function in p5.8p . tt® in p5.9h and sum up the equations. 


—[ C(A' 
12 2 dtj^ ^ 


'p(u)) : (p(u)A'^)-v/adzidz 2 + 1(1 zsTr^dzo^A'^ : p{u)^/adzldz^ 

X + 2 Juj \J~l/2 J 


1^ 
2 dt 


—A'^ ; 7 — z^^J^ —A'^ : p(u) ) Tr^y/adz 
A+2 A+2 ' 


f / "X ^ /* ^u. C 

+ J '/adz = J {V++V-) ■ —/adzidz 2 ^ J V//ads. 

The anti-symmetric terms cancel out as before. We integrate the equation over time and use 
the initial conditions to obtain 


1 1 
12 2 


/ C(A'^p(u)):(p(u)A^)Vadzidz 2 + ^ [ (^ + ^ {/)'^/adz 

Jui 2 \ A + 2/ 

, f f 

Jo Jn V 9z3 j Jo Jq 


a- - A^:-^//adzdT 

n X + 2 dt 


(5.13) 


f f {V++V-) ■ ^/adsdr T [ [ VZ/adsdr. 
10 Jlu dt Jo 
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Next we compute the elastic energy 

[ C (Q7°Q^) : (Q7°Q^) V^dz 
Jn 


= / A(tr (Q7°Q'^))^ + 2 Q’^Q 7 °Q^Q : 

Jn 



= / A 


— —A'^ : 7 — zs^^ —A'^ : p(u) + — 7 r° ) ^dz 

n VA+2 A+2 A+2 ' 



+ 2 ( A ^(7 - Z 3 p(u)) A'" : (7 - 2 ; 3 p(u)) + 


-TT — 


A+2 A+2 


A'" : (7 - 2 : 3 p(u)) ]y/ddz 


= 1A(2. " 



n \A+2A+2 


vr^A*^ : 7 - 2 


+ 


\ A + 2/ 


(A'^:7)' + (^3)' 


A+2A+2 

V A + 2/ 


tt^zsA'' : p(u) 


(A'^ : piu)f + 


(l^) 


+ J 2 |^A'= 7 A‘' : 7 + (2:3) A‘'p(u)A'= ; p(u) 

_ 2 , ^ ' 7 r°A" : 7 + 2 , ^ ^ 7 r° 23 A" : p(u) 


A+2A+2 


+ 


V A + 2/ 


^0\2 _L 

TT ) + 


A 


A+2A+2 

2 


(A^ ; 7)2 + (23)2 f j (A^ : p(u))2^Va<^2 


A(2)2 2A2 


/nVV(A + 2)2 (A+ 2)2 


A + 2y 

(A'^ : 7)2 + 2 A'^ 7 A'^ : 7 


+ (^ 3 )' 


+ A 


+ 


2A2 


(( A(2)2 

l^l^(A + 2)2 ' (A + 2) 


(A'=:p(u))2 + 2A'^p(u)A‘=:p(u) 


a 


+ 2 - 


a 


(A+ 2)2 (A+ 2)2 


+ 



2Aa- 


— 4a- 


( 7 r °)2 j ^/adz 

A 


+ I ~A2- 


(A + 2)2 (A+ 2)2^ 

a 2 , a A 

+ 4. 


A+2A+2 A+2A+2 


vr^A'^ : 7 

) vr^zsA"^ ; p(u) ) y/adz 


C{A^'r) : + (z 3 ) 2 (:(A'=p(u)) : p(u)A‘' + ) Vadz. 


Insertion of the above equality and (15.1311 in ()5.12p . yields 


1 1 


^ = / C(AV(u)):(p(u)A‘=)v^d2;id2:2 + 7 / /3 + 


12 2 


+ 


a 



10 Jn 


dTT^ 

dz?. 


y/adzdr + / 

Jo Jn 


2 Jn V ' a + 2 

^ 2 AC ^7 0 J 

a- - A -.—'K yjadzdr 


(tt^) y/adz 


A + 2 dt 


C(A'=7):7A‘' + {zzfC{A^p{\i))-.p{\i)A‘' + ) ^dz 

F) 0 \ 2 

y/adzdr 


b [ {TryVddz - [ [ ( 
^ Jn Jo Jn V 


5z3 
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fl 


2 AC ^ 

-A'": — 

fl A + 2 dt 


'K^^/adzdT ~ 2 j V^dz. 


From 


for q independent of transversal variable we obtain 


- f 

dt Jfi 


f f /3 -h ^ vr° + : 7^ q\fadz = 0, g G H^{uj). 

\\ A + 2/ A + 2 J 


This implies that 


a- 


A + 2 dt 

Inserting (I5.14p into A yields 


A':A = -|/,+ 


a 


4 vrOrfzg. 


A+ 2/ dtJ_i/2 


(5.14) 


A = - /3 + 


a 



d 


»l/2 


A + 2/ 7o Jn dt \ 7-1/2 


TT^dzs ) TT^y/adzdr — - (C(A'^7) : 7A‘^) y/adz 


= - /3 + 


a 


1 


= - /3 + 


A + 27 2 7o 

,2 


dt 


Lie 


TT^dzs y/adzidZ2dT — 


C(A'= 7 ) : 7 An y/dd. 


a 


A + 2 


MZ 


TT^dzs J y/adzidz2dT — - (C(A^7) : 7A‘^) y/adz 


(5.15) 

where in the last equation we have used that vr^|t=o = 0. Since A > 0 by definition and since 
the right hand side is nonpositive we conclude that A = 0. Positivity of C implies 7 = 0 and 
thus the strain 7 ^ is fully determined by u and 


7° = 


-Z3p{u) 


0 

0 


0 0 Afi-° + "3x^A-p(u) 


Moreover, from (I5.14p we obtain that = 0 thus the poroelastic flexural shell 

model is given by 


Z f C(A‘'p(u)) : p(v)A''v/adzid 2;2 + 


1/2 


X 2 JU! —1/2 

= / {V+ + V-) ■ Vy/adzidz2, v G Vir(w). 

J U1 

L ('’+ jZ) ■ L ^ + 

Vqy/ads, q G 


z^TT^dz^A'^ : p{v)y/adzidz2 


(5.16) 


f dir^ dq 
Iq dz3 dz3 


y/adz 


(5.17) 


= + 


We now have A(e)(t) —)• 0 for every t G [0,T]. Since A(e) : [0,T] —)• M is an equicontinuous 
family, we conclude strong convergences of the strain tensor and the pressure 


1 


7 ^(u(e)) ^ 7 ° 


7r(e) —7r° 
dTT{s) dir^ 


dz3 


dz3 


strongly in (^([O,T];+"^(12; 
strongly in C'([0, T];+^(12)), 
strongly in L^(0, T; L^(12)). 


(5.18) 
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Step 4 (Strong convergence for displacements). We prove the strong convergence in two 
steps. In the first step we use the last part of Theorem 1131 and prove pointwise convergence 
of u(e). Due to equicontinuity it then implies the uniform convergence, i.e., we obtain 

u(e) —^ u strongly in (^([OjT]; 

Lemma 17. 

p(u(e)) —^ p(u) strongly in (^([O,T]; 

Proof. From [71 Theorem 5.2-1] for a,/I = 1,2 we have the estimate (for a.e. t G [0, T]) 

1 d 

l|-^7a/j(u(e)) + Pa/3(u(e))||^^-l(t^) 

3 (5.19) 

< llT^3(u(e))||L2(o) + ^II^i(^)IIl2(o) + £||'*^2(e)||L2(o) + e||^^3(£)||_H'i(o)), 

i=l 

for C independent of e. By the strong convergence (15.181) scaled transformed symmetrized 
gradient i 7 ^(u(e)) is bounded in (7(10, Tj; L^(D; M^^^)) and by the a priori estimates from 
Corollary [12] u(e) is bounded in (^([O, T]; i7^(D; M^)). Therefore the right hand side in (|5.19p 
tends to zero uniformly with respect to t G [0,T]. Thus we obtain 

-7^7a/3(u(e))+ p«/3(u(e)) ^ 0 strongly in (7(10, Tj; i7'^(D)), a,/3 = l,2. (5.20) 

e 0 Z 3 

From (I5.18P we have that 

-^7"(u(e)) ^ ^ 7 ° strongly in C{[0,T]; 
e dz 3 dz 3 

Using this convergence in (I5.20p we obtain that functions p(u(e)) converge strongly in 

(7([0,r];F-i(D;M2x2))_ □ 

Theorem 18. 

u(e) —)■ u strongly in (7([0, Tj; M^)). 

Proof. From Lemma fT71 and Remark [TT] we have pointwise convergences 

p(u(e))(t) —^ p(u)(t) strongly in 77“^(D; 

u(e)(t) ^ u{t) weakly in i7^(D;]R^), 

for every t G [0, T]. Thus the last part of the Theorem [13] (taken from [?]) implies 

u(e)(t) ^ u(t) strongly in i7^(D;M^), 
for all t G [0, T]. This implies that the function tc(e) : [0,T] —)• M given by 

w{e){t) = ||u(e)(t) - u(t)||^^i(n;R3) 

converges pointwisely to zero. From the uniform estimate of u(e) in 17^(0, T; iL^(D; M^)) we 
obtain equicontinuity of the family (r(;(e))£>o. This implies the uniform convergence of w{e). 
This implies the statement of the theorem. □ 
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Step 5 (Strong convergence for the stress tensor). As a consequence of the convergence 
of the term i 7 ^(u(e)) from (j5.18l) we obtain the convergence of the scaled stress tensor. 
Proof of Corollary [3 We compute a in the local basis given by Q using (15. 5p and (|5.6I) 
and that 7 = 0. We obtain 

Q^aQ = Atr (Q 7 °Q^)Q^Q + 2 Q^Q 7 °Q^Q - avr^Q^Q 


aX 


2A 


-TT - 2 : 3 - 


+ 2 


A + 2 A + 2 

— 2 ; 3 A'^P(u)A^ 


2 a 

A + 2 


0 


TT - 23 



2A 

A + 2 


A^:p(u) 


A^ 0 
0 0 


A+2 


(^(A'':p(u))I + 2A''p(u 


^vr0 + 23 ^A+p(u)) 


223A'=p(u)A‘= 0 
0 0 


0 

0 


□ 


6 Appendix 


6.1 Properties of the metric tensor, the curvature tensor and the third 
fundamental form 


Some symmetry properties of geometric coefficients are listed in the following lemma. For 
the proof see [ 6 ] . 

Lemma 19. The following symmetries hold (a,l3,K G {1,2}) 

aag = a Pa, Kp = bpa, ^a/3 = 


The change from the basis to basis is done using 


Moreover, one has 


2 2 

K,= l K =1 

2 2 

biy = b}\a, = 

K = 1 K.= l 


6.2 Computation of 0(v) 


Lemma 20. For v G Vf{(^) and ©(v) defined in the Step 2 of the convergence proof (Sec¬ 
tion 0) one has 


0(v) = -Z 3 


P(v) 


2 

J^biid^vs + Vrbl) 

K, = l 


2 

K, = l 


2 

J^biid^vs + Vrbl) 

K = 1 
2 

J2b2{dKV3 + Vrbl) 
K.= l 

0 
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Proof. Let us first denote 


= = E 


K,= l 


dii^divs + 2v^,b1) i {d2{^diV3 + 2v^bf) + 81(^821 

{82{^8iV3 + 2vK.b1) + 8i{^82V3 + 2Vf,b'^)) 82{^82V3 + 2vKb'^) 


Then using (|4.3p and (|4.4p we obtain 


Q(v) ^ 23 [ - s + '^{d^V 3 + 2 ^ ■ 


+ E' 

K =1 


■pK 

J- 11 

PK 

21 


-pK 

12 

pK 

22 


-bf 

— b^ 
“2 


K,= l 


T=1 

_ -b1 

-b1 

0 


bill 

b1\i 

h\b1 ' 

2 

■ blb^i 

b1b^2 

0 

b1\2 

b1\2 

blbl 

+ V3^ 

b^bnl 

&2^k2 

0 

blbl 

hlbl 

0 

K=l 

0 

0 

0 


= -^3 


E 

, K,= l 


^8(v 3 + 28iVi^bf 
I 8182 V 3 + 82 V^bf + 8 iv^b'^ 
0 


8182 V 3 + 82 V^bf + 8 iVnb^ 
+ 282 V 


0 


2 

28 ib 1 

82h1 + 8ib1 

0 ■ 

2 

r 

^ 11 

pK 

3 12 

0 ■ 


82h1 + 8ib1 

282b1 

0 

- y] 8i^v3 

pK 

^ 21 

pK 

3 22 

0 

K=l 

0 

0 

0 

K =1 

0 

0 

0 


2 


pr pr 

3 11 -'-12 

0 ■ 

2 

■ bill 

& 2 I 1 

0 ■ 

2E 

VKbl 

pr pr 

3 21 -*-22 

0 

-iZ^K 

b1\2 

& 2 I 2 

0 

K,,T=1 


0 0 

0 _ 

K,= l 

0 

0 

0 _ 

2 

‘ blbi^i b1bi^2 

0 ■ 

2 

2 




^^3 


E 

«;=! 


^ 2 ^Kl 

0 


^ 2 ^k 2 

0 




«:=! 


r=l 


0 

0 

"2 


-6^ 

— h^ 
"2 


Expressing 9a using 6 ^|q and collecting all terms with we obtain 


0(v) = - Z 3 \ 


/ 4 , 

l 8 fv 3 + 28 ivJ )1 

i 9 i 92 U 3 + 82VJ)1 + 8ivJ)1 0 ‘ 


\8182v3 + 82V J)1 + 8ivJ)1 

\8lv3 + 282V 0 

^ K =1 

0 

- 1 

0 

0 


+E' 

«:=! 


&?li- 2 E'=ir?A 
E"=ir 2 ^ 6 I + 62 li-E?=i 


b^\2-El=i^^rbl-El=i^ 


pK LT 

^ lr"2 


65 I 2 - 2 EL 1 T 


K ip- 

1t"2 


^ Z-yT=l - 
0 


2.^2 


- X] 

2 

+ ^(9^i;3 + ^Ur 6 ;; 


r 

-*-11 

pK 

3 12 

1 

0 

2 


b1b^2 

1- 

0 

pK 
-*- 21 

pK 

-*- 22 

0 


blbnl 

b2bK2 

0 

1- 

0 

0 

0 _ 

K = 1 

0 

0 

0 _ 


k ;=1 


T = 1 


0 

0 

hi 


0 

0 

hK 

"2 


b1 

UK 

"2 

0 


This implies the statement of the lemma. 

A part of the proof can be also find in (Tj Step 4 in Section 6.2], 


□ 


6.3 Cylindrical surface 

Let ujl = {~L/2,L/2) x (0, d), where d G (0,27r) and let {z,6) denotes the generic point in 
u. Let i? > 0. We define the cylindrical shell by the parametrization 

ip : ujl ^ (p{z,9) = {RcosO, RsinO, z)"^. 
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For d = 2 tt the surface is the full cylinder. Then the extended covariant basis of the shell 
S = (p{uj) is given by 


ai{z,9) 

a2{z,e) 

a3{z,e) 


d^(p{z,e) = ( 0 , 0 , 1 )'^, 
dg(p{z,9) = {—Rsin9,Rcos9,0)^ 
ai{z,9) X a2{z,9) 

I / D\ / D\\ ^ COSu, 

\ai(z,d) X a2{z,9)\ 


sm9,0f. 


The contravariant basis is biorthogonal and is defined by 

a\z,9) = (0,0,1)^, 

a^{z,9) = (-^sin 6 », ^cos 6 l, 0 )^, 

a^{z,9) = {—cos9, — sm9,0)^. 


The covariant Ac = (oq,/?) and contravariant A'^ = (a“^) metric tensors are respectively given 
by 


■ 1 

0 


■ 1 

0 

0 

R^ 

, = 

0 

1 

w \ 


Ac = 

and the area element is now ^/adS = Vdet A^dS = RdS. The covariant and mixed compo¬ 
nents of the curvature tensor are now given by 

fell = &12 = ^21 = 0, b22 = R, b\ = bl = bl = 0, bl = ^. 

A simple calculation shows 


= a,/3, cr G { 1 , 2 }, 

fe^U = dab’p — = dab^p = Q, a,/3, u G {1, 2}. 

Now the displacement vector v in canonical coordinates is rewritten in the local basis 
V = Qv = vio^ + V 20 ? + v^a^. Note that contravariant basis is different than the usual basis 
associated with the cylindrical coordinates. One has 


vi = Vz, V2 = RV0, V3 = -Vr. 


Similarly, P± = Q ’^P± = {V±)iai + {V±) 2 a 2 + ('P±)3a3- Thus 


Thus 


{P±)l = {P±),, {P±)2 = ^{R±)e^ {'P±)3 = -i'P±)r- 


V±-v = {V±)iVi + {V±)2V2 + iP±)3V3 = {V±)^Vz + {V±)gVe + {V±)^ 


Inserting the geometry coefficients into the strains 7 and q we obtain 


7(v) = 


<91^^! - ELi - fellll3 

\{dlV2 + d2Vi) - Ek =1 '^21'^K - b2lV3 
dzVz 


\{diV2 + d2Vi) - Y?k=1 bl2^^K - fel2?^3 

d2V2 - ELi ^22Vk - b22V3 


p(v) = 


L ^{RdzVg + dgV^) 

dnV3 di2V3 
d2lV3 d22V3 
2 

+E 


^{RdzVg + dgVz) 
RdgVg + RVr 


K=1 


bidiv,^ + bldiVi^ + b1\iVi^ - fe^fe^ius b^div^ + feJ52U«, + b1\2V^ - b^b^2V3 
bid2V^ + b^div^, + fealiu^ - b^b^iV3 b'^d2Vf, + b'^d2V^ + b^\2V^, - b'^b^2V3 


-dzzVr -dzOVr + dzVg 

-dezVr + dzVe -dggVr + 2dgvg + Vr 


29 


















As example we write the model on the space 

= {{Vz,Ve,Vr) G X H^{uJl) X H‘^{uJl) : {Vz,Ve,Vr)\B=Q,d = 0 , 

1 ( 6 . 21 ) 

deVr\e=o,d = 0 , = -{Rd^ve + dev^) = Rdeve + Rvr = 0 } 

which includes clamping boundary conditions only on two generatrices of the portion of the 
cylinder. For cylindrical shell fully clamped one has = {0}, and the shell behaves as 

the generalized membrane shell, see [TJ Section 5.8]. From the condition of inextensibility in 
(| 6 . 2 ip we obtain 

Vz{ 0 ), d^V 0 = -^dev^, Vr = -dgve. 

Therefore (using notation ' = dg) 

v^{d), ve{z,e) =-^v'^{0)+ wb{ 0), Vr{z,9) = -w'g{9). 

Smoothness and boundary conditions for Vz,vg,Vr imply 


Vz£H^iO,d), wg€H^{0,d). 


Thus 


p(v) = 


0 

-W: + VzY -Wz + v.r + « + weY 
Now we insert this into the model given by ()3.2ip ~ (l,'1.23p written in dimensional form 

f ( 2/rA 1 

12 



0 J-LI2 


+ (J 4 + ^zY{v" + v^Y + + ^z:Y' - i^e + ^eY){^iv" + v^Y' - i^g + wgY)^ 


fd /■V2 rtj'i „ 1 / Z \ 

/ / rp^dr^ ^ '^dY)) Rdzd9 

0 — Lj A •J —t/z 


+ 


2pa 

A + 2p Jo J-L/2 J-l/2 

d fLl2 


0 J-L/2 


nX + YPl%)^z + («)^ + YPl%){--v', + wg) 


+ + (ip-^),)(|< - w'g))Rdzd9, V, G H^{0,d),wg G Hl{0,d). 


d rd ,L/2 

dt J-e /2 Jo J-L /2 


I3g + 


a 


A -|- 2 /i 


p^qRdzd9dr 


2ap /'^/2 rd rLl2 -y 


+ 


A + 2 // 7_£/2 Jo J-L/2 R^ ^ R 

k dp^ dq 


^ (“(^^*('“2 + '^zY' + dtiujg + ujgY^ rqRdzd9dr 


k r r 

V J-i/2 Jo J-L/2 


dr dr 


Rdzd9dr = 


d rL/2 



0 J-L /2 


^L(g(-x) - q{7;))Rdzd9 


in V'{0,T), qGH\{-il2,il2y,L\iVL)), 

= 0 at t = 0. 

After integration over 2 the first equation separates and we obtain the following problem: 
hnd {p^ ,Uz,ujg} G C{[0,T]-, x HQ{0,d) x HQ{0,d)), drP^ G L^((0,T) x satisfying 


^ Rdzd9 


30 















the system 


rd ^ 

Jo ^ 


12 In i?3 




X + H 

A + 2/i 12R‘^ 


« + u,nv'^ + V,)" + 4/xL« + u,)'{v'l + V,)' dO 


2 //« 1 


A + 2/i Jo \J-L/2 J-e/2 

rd / fLl2 


/ / rp^drzdz {v” + Vz)"d9 

J-L /2 J-e /2 I 


rd / /-i/Z /■i/2 

[R {{Rt^ + {Vl%)dzv, - / {{■Pt% + {VZ%)zdzv', 

Jo V J-L/2 J-L/2 

r-L/2 \ 

/ {{Rt^ + {Rl\)zdzv'i de, V, S H^iO, d), 

J-L/2 J 


rL/2 
^-L/2 

t r L ^ 
12 Jo R^ ^ 


+ 


I ( Ll, Ll, ’■*’"**) <“»+ 

/■i/2 /■i/2 \ 

/ + / ((^^')r- + {Rl%)dzw'g Rd9, we G ^^(O, d), 

'-L/2 J-L/2 / 


dt J-e /2 Jo J-L /2 


Pg + 


a 


A J- 2// 


p^qRdzd9dr 


2ap 1 


/. 


+ 


A + 2/1 J _£/2 Jo J-L /2 


A: r r 

V J-e/2 Jo J-. 


RJ rqRdzd9dr 

d r-L/2 


Rdzd9dr = 


yL{q{--) - q{-))Rdzd9 


-e /2 Jo J-L /2 dr dr jq j_l /2 

in P'(0,r), gG JjH(-^/2,^/2);L2(a;L)), 

= 0 at t = 0. 

The terms in the shell equation appear in the classical model of linear model of cylindrical 
shells, see e.g. [E] and d]. 
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